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EXPERIMENTAL AND ANALYTICAL STUDY OF THE TRANSIENT SOLIDIFICATION 
OF A WARM LIQUID FLOWING OVER A CHILLED FLAT PLATE 
by Joseph M. Savino and Robert Siegel 
Lewis Research Center 

SUMMARY 

Analytical predictions and experimental data are provided for the transient growth of 
a frozen layer that forms when a warm liquid flows over one side of a chilled flat plate. 
The plate is cooled to a temperature below the freezing point of the liquid by a coolant 
flowing along the other side of the plate. This could be encountered, for example, in a 
heat exchanger where water is being cooled by a cryogenic fluid. Three analytical pro- 
cedures are employed and compared for accuracy and convenience in application. One is 
an iterative technique that yielded three closed-form solutions, each successive solution 
being a more accurate higher order approximation. Numerical results from the third 
approximation are presented graphically in sufficient detail so that rapid estimates of 
frozen layer growth can be made. In the experiment, water flowed over one side of a thin 
Inconel plate while the other side of the plate was cooled by flowing chilled alcohol. Tran- 
sient ice growth was initiated by a sudden increase of the alcohol flow. The ice thickness 
and the plate surface temperatures were measured throughout the transient growth period 
and compared favorably with theoretically predicted values . 


INTRODUCTION 

This report is concerned with the transient solidification of a warm flowing liquid that 
is in contact with a chilled plane wall. The wall is cooled below the freezing point of the 
warm liquid by a coolant flowing on the opposite side. As the frozen layer grows on the 
surface , the warm liquid is continually supplying convective heat to the liquid - frozen 
layer interface. This convective heat combines with the latent heat released at the inter- 
face and is transferred through the frozen layer and wall into the coolant. This type of 
solidification problem arises in important and familiar applications, such as continuous 
casting of metals and freezing of rivers. The application that motivated the present study 



is concerned with a number of advanced propulsion devices that utilize liquid - liquid heat 
exchangers in which the coolant is at a temperature below the freezing point of the warm 
liquid. One example is a water -to -liquid -hydrogen heat exchanger for possible applica- 
tion in a nuclear rocket engine. If the temperature or flow rate of the warm liquid is low 
enough, a frozen layer will form on the exchanger walls. Such an exchanger can operate 
satisfactorily if the frozen layer is not thick enough to severely restrict the flow passages. 
The prevention of serious blockage due to excessive solidification requires an understand- 
ing of the transient formation of the solid layer on the exchanger walls. 

A review of the literature was made to uncover what information was available on so- 
lidification or melting problems with a convective boundary condition at the moving inter- 
face. Little analytical work was found for this type of boundary condition. For other types 
of boundary conditions, however, there were numerous mathematical studies as evidenced 
for example by references 1 to 5. Relatively few experimental studies were found, and of 
these none were for forced flow of the liquid over the moving interface. Some of the ana- 
lytical studies are now reviewed since they form a foundation for the present analysis. 

There are two reasons motivating the considerable number of mathematical studies. 
The first is the practical importance of moving interface problems in applications such 
as thermal storage devices for space vehicles, the formation of polar ice, and the freez- 
ing and thawing of soil. The second is that the problems offer a challenge for mathe- 
matical study because the moving boundary introduces a nonlinearity into the mathema- 
tics. The nonlinearity considerably increases the difficulty of finding solutions. As a 
result, exact solutions for transient frozen layer thicknesses and the accompanying tem- 
perature distributions and heat flows are usually difficult to obtain. Exact analytical so- 
lutions have only been found for simple configurations with certain types of boundary con- 
ditions. This has led to the situation where, in some instances, physically unrealistic 
boundary conditions have been imposed to acquire a mathematical solution. Some exact 
solutions for conditions of practical importance are summarized by Carslaw and Jaeger 
(ref. 6). One of these, for example, is the solution by Neumann for freezing (or melting) 
of a stationary semi-infinite medium initially at constant temperature, such as freezing 
a lake or melting a thick metal piece without removing the molten material. The semi- 
infinite medium is being cooled (or heated) by maintaining the boundary at a constant tem- 
perature. 

When an exact analytical solution cannot be found, numerical or approximate methods 
must be employed. Murray and Landis (ref. 7) devised a numerical finite difference tech- 
nique and demonstrated it for the following sample problems: 

(1) A semi-infinite medium initially at the fusion temperature is subject to a step re- 

duction in temperature at its boundary 

(2) A slab of molten material of finite thickness initially above the fusion tempera- 

ture is solidified by having its boundaries suddenly become subcooled. 
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Cochran (ref. 8) developed a lumped -parameter method wherein all the heat capacity 
is lumped at the center of the body and applied the method to a slab, a cylinder, and a 
sphere. A heat-balance integral approach was demonstrated on a variety of freezing and 
melting problems by Goodman (ref. 9). The temperature distribution in the material was 
approximated by a second- or third -order polynomial to evaluate an integrated form of the 
conduction equation. Adams (ref. 10) employed an iterative technique to find the frozen 
layer thickness and heat transfer in plane-, cylindrical-, and spherical-shaped metal 
castings with and without accounting for the conduction into the mold. This method does 
not involve any approximations in principle and provides successively iterated solutions 
that should approach an exact solution. In some instances, only the first iteration was 
conveniently carried out in analytical form. Cullom, et al. (ref. 11) employed the Crank - 
Nicolson finite difference technique to solve the problem of solidification of a stationary 
liquid at the fusion temperature onto a round tube cooled from within. 

As mentioned, experimental results relevant to the present study are quite limited. 
Thomas and Westwater (ref. 12) studied freezing and melting in a thermal-conductivity- 
type cell. The cell was heated at the top, cooled at the bottom, and insulated on the sides 
to provide a unidirectional heat flow by conduction only through the freezing liquid. The 
primary test substance was n-octadecane and, for all conditions, the solid-liquid inter- 
face was seen to be populated with microscopic peaks and valleys. The measured veloci- 
ties of the solidification front were as much as 100 percent higher than theoretically pre- 
dicted. Observations of melting and freezing of p -xylene and water revealed that these 
substances provided smooth interfaces in contrast to the crystalline interface found for 
n-octadecane. 

Tests of freezing and melting of water in absence of forced flow were conducted by 
Boger (ref. 13). Water was contained in a chamber 2 inches high bounded on top and bot- 
tom by copper blocks. The temperatures of the blocks could be regulated to obtain freez- 
ing or melting at either the upper or lower block. For some tests, the temperature range 
of the water within the chamber was entirely between the freezing point and the density in- 
version temperature (4° C). With these conditions, tests were conducted with no buoy- 
ancy forces present, and the freezing front velocities were found to agree quite well with 
values predicted by using the numerical technique in reference 7. When free convection 
was present, it was necessary to modify the theory by introducing an effective thermal 
conductivity to account for the convective transport. 

Some experiments for melting in the presence of free convection are given in refer- 
ences 14 and 15. Reference 14 studies the melting of ice cubes by free convection in a 
water bath. In reference 15, the lower end of a vertical cylinder of ice was melted from 
below in a one -dimensional fashion by free convection from heated water. Comparison 
with theory indicated that the ice front melted more slowly than analytically predicted. 

Solidification inside a tube and channel in the presence of laminar forced convection 
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was studied both experimentally and analytically by Zerkle (ref. 16). The wall of a cir- 
cular tube and a pair of opposite walls of a rectangular channel were cooled so that ice 
layers formed. Zerkle was primarily interesting in developing a model for predicting the 
steady -state heat transfer and pressure drop inside channels whose walls were coated 
with ice. 

In analytical models, the interface is usually assumed to be smooth and at the equi- 
librium freezing temperature. The details of crystal structure and undercooling at the 
solid-liquid interface have been studied in a number of experiments (refs. 17 to 21). For 
ice, if the propagation is into supercooled liquid where the heat of fusion is being dis- 
charged into the liquid, the interface is covered by a lattice of protruding crystals. For 
the present experiments, the liquid is above the freezing temperature, and the heat of 
fusion is being withdrawn by conduction through the ice layer. For these conditions, a 
smooth ice layer would generally be expected. Some comments concerning the interface 
will be made later in the discussion of the present experiment. 

The present investigation was undertaken because of the lack of both analytical and 
experimental information on the transient solidification in a flowing medium. The objec- 
tives were to devise an analytical method for predicting the transient thickness of the 
solidified layer and accompanying temperature variations when a flowing fluid freezes on 
a cold plate and to compare the analytical results with experimental findings. 

The model chosen for the analytical study is a locally one -dimensional solidification 
of a liquid that flows over one side of a thin flat wall. The other side of the wall is being 
convectively cooled by a fluid that is at a temperature below the freezing point of the liq- 
uid. The wall on which the liquid solidifies is assumed to be sufficiently thin so that the 
heat flow needed to subcool the wall during the transient is negligible compared with the 
heat flow through the wall. The heat flow through the wall is comprised of the latent heat 
of fusion, the subcooling heat capacity of the frozen layer, and the convection from the 
warm liquid at the solidified interface. 

Three analytical methods were employed. For two of these, the transient heat- 
conduction equation that governs the energy balance within the solidified layer was inte- 
grated in a general way to provide an expression for the time required to form a speci- 
fied instantaneous frozen layer thickness. The general integrated equation is evaluated 
for the present problem by using a method of successive analytically iterated approxima- 
tions and also by using an approximate temperature profile of the type proposed by 
Goodman (ref. 9). The iterative technique is applied in this study in a different manner 
than the way it was applied by Adams (ref. 10). By the present method, higher order 
solutions can be found, and an analytical expression very close to an exact solution is ob- 
tained. The third analytical method is one proposed by Goodman (ref. 9) in which a heat 
balance at the liquid - frozen-layer interface is used. 

The methods from reference 9 yielded two solutions that had to be integrated numer- 
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ically. The iterative method yielded three closed-form analytical expressions by means 
of three successively improved iterations. All the solutions are compared with regard to 
accuracy and ease of application. One of the closed-form solutions is sufficiently simple 
that it can be readily evaluated, and it yields quite accurate solidification times for all 
ranges of the parameters. Another closed-form solution, of a more complex form, rep- 
resents for all practical purposes the exact solution and is used to evaluate the final re- 
sults of the analysis that are presented in graphical form. 

In the present experiment, ice was formed on a thin Inconel plate. Water at a bulk 
temperature above the freezing point passed in forced flow over one side of the plate while 
the other side was cooled by flowing chilled alcohol. Transients were initiated by sud- 
denly starting the alcohol flow. The transient variations of ice thickness and plate surface 
temperatures were measured and are compared with theoretical predictions. 


ANALYSIS OF SOLIDIFYING LAYER 

The analysis deals with a one -dimensional configuration, as shown in figure 1. This 
represents an idealization of what can be achieved experimentally. Although the present 
experiment is for an ice -water system, the analysis could apply to any fluid that behaves 
within a reasonable approximation according to the conditions assumed (e.g. , constant 
properties, solid-liquid interface at the equilibrium freezing temperature). 

As a prelude to the analysis, it is instructive to trace the solidification process qual- 
itatively. Before the start of any transient, a warm liquid at a fixed temperature T^ 
flows over one side of a thin plane wall. The warm liquid provides a convective heat- 
transfer coefficient h^ that is assumed to remain constant throughout the entire tran- 
sient period before and during solidification. A coolant at a fixed temperature T c flows 
over the opposite side of the plate where it provides a constant convective heat-transfer 
coefficient h c . The resulting steady heat conduction through the plate establishes a 
linear temperature profile within it. If no coolant is flowing, the plate is at the tempera- 
ture T^ of the warm liquid. 

A transient solidification process can then be initiated in a number of ways, such as 
reducing the temperature or flow rate of the warm liquid, introducing the coolant if it had 
not been flowing, or lowering T of the flowing coolant. These changes cause the wall 
to further cool until the surface on the side of the plate exposed to the warm liquid 
reaches the freezing temperature. At this time, since the specific heat of the wall is 
neglected as an analytical simplification, a linear temperature distribution exists in the 
plate, as shown by the dotted temperature profile in figure 1. Solidification is assumed 
to begin as soon as the plate surface reaches the freezing point; that is, any subcooling 
necessary to initiate nucleation is neglected. 
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As the solidified layer on the wall increases in thickness, heat is extracted from the 
boundary layer of the warm liquid and is transferred by convection to the liquid-solid in- 
terface. There, it combines with the latent heat released by the solidification process and 
and is conducted through both the frozen layer and cold wall to the wall surface in contact 
with the coolant, where it is convectively removed. In general, an additional amount of 
heat removal is needed to subcool the frozen layer, and this is also transferred to the 
coolant (the heat removal necessary to subcool the wall is neglected in the analysis in 
comparison with the other heat flows). The solid layer continues to grow until it achieves 

a steady-state thickness X . 

s 

In the analysis and experiment, the time t required to grow any given layer thick- 
ness X(t) < X_ is measured from the moment solidification begins, which is taken as 
t = 0. The properties of the wall and the solid layer are assumed to be independent of 
temperature. 


Steady-State Thickness of Frozen Layer 

In many solidification problems, the frozen layer never approaches a steady-state 
thickness; rather it continues to propagate with time into the liquid. When, however, the 
liquid flowing over the frozen layer is at a temperature above the freezing point, the 
frozen layer will achieve a steady-state thickness. This thickness will be used as a ref- 
erence length in the later analyses of growing layers and hence will be separately derived 
here. 

It is assumed that the liquid-solid interface at the boundary of the frozen layer is at 
the freezing temperature Tj. Taking heat flow as positive in the positive x-direction, 
the convected flux from the liquid to the liquid-solid interface is 

-q = hj (Tj - Tj) s constant 


Using the overall heat-transfer resistance between the liquid-solid interface and the cool- 
ant gives the relation 


-q = h, (T ; - T f ) 



k +k w\ 


This is solved for the steady -state thickness 


6 



( 1 ) 




For X = 0, equation (1) gives the relation between variables required to just avoid 
s 

freezing. For example, solving for gives 


z| x s =° 



( 2 ) 


In order to prevent freezing, the liquid temperature must be equal to or greater than 
the value given by equation (2). 


General Equations for Frozen Layer Growth and Temperature Distribution 

The heat flow within the frozen layer is governed by the transient heat- conduction 
equation; since constant properties are assumed, this has the form 


ax 2 P 3T 


(3) 


When equation (3) is integrated from any position x within the layer to the solid-liquid 
interface X, there results 


k 




(4) 


At the interface, the heat conducted into the solidified layer is equal to that supplied by 
the latent heat of fusion and the convection from the flowing liquid; that is, 

= P L^+ h 7 (T 7 - T f ) (5) 

dr L L 1 



Equation (5) is substituted into equation (4) to yield 
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-k®* 

dx 


X 


( 6 ) 


= - pL f -M T ; -w 

T <X 


0T 

0T 


dx 


The term on the left of equation (6) is the heat flux crossing any position x at any time t. 
The first two terms on the right are the latent and convective heat fluxes originating at 
the solid-liquid interface, while the last term is the heat being removed by subcooling the 
portion of the solidified layer between x and X. The integration of equation (6) from the 
wall (x = 0) to any position x results in an expression for the instantaneous temperature 
distribution 


T(x, t) - T, = £LdX x + HL(t . T w -fi f K f X 
2 k dr k 1 k J 0 J x 


0T 

0T 


dx dx 


( 7 ) 


To eliminate the variable Tg, which is a function of time, equation (6) is written at 
x = 0 to give the heat flow at the wall, and this is equated to the heat flow through the wall 
and into the coolant to provide the relation 


q(x = 0) = -pL ^ - h z (T z - T f ) + pc p dx 


dr 


•'0 


0 T 


T 2 ~ T c 
1 a 

h c +k w 


( 8 ) 


Then equation (8) is solved for Tg and is substituted in equation (7) to give 


T(x,t) = T + (J- + -L. 

V h c k w / 


P L ® + \ < T Z - T f> - P', 


f X 0T 
JO ^ 


dx 


+ ^ — x + — (T 7 - T f )x 


k dr 


pep /.x rX 

k Jo Jx 


0T 

0T 


dx dx (9) 


Equation (9) is placed in dimensionless form by letting 
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h ; ( T l - T f) r 


T - Tr 


T' = 


T c- T f 


and using equation (1) in the form (X g h z /k)[(T^ - Tj)/(T f - T c )J = R/(l + R) to give 


i - T* = 1 + **'(<*: ♦ ii ♦ r h; dx' + r f x ' ax' dxi 

1 + R Wo 3T ' Jo Jx' dT ’ / 


1 + Rx'i 

f dX' 

1 + R 

Idr* 


( 10 ) 


Because of the nondimensionalization, it is noted that two parameters, 

S = Cp(T f - T c )/L and R = (X g /k)y/ [(l/h c ) + (a/k w ) , have appeared. The subcooling 
parameter S is the ratio to the latent heat L of the maximum internal energy that could 
be removed if a unit mass of the solidified material were subcooled from the freezing 
temperature to the lowest temperature in the system T . This is not, however, a mea- 
sure of the subcooling that actually occurs during the freezing process, since the entire 
solidified material does not reach this low temperature. The parameter R is a ratio of 
heat flow resistances. The X /k is the resistance offered by the steady-state thickness 
of the frozen layer, while (l/h c ) + ( a A w ) is the combined resistance offered by the wall 
and the convective coefficient on the coolant side of the wall. 

By applying the rules for differentiating under an integral, the integral terms of 
equation (10) can be transformed so that the derivative is taken outside the integral signs 
and equation (10) becomes 


T’ 


1 + Rx'/dX' + 
1 + R \dT’ 


RS 

1 + R 



T’ dx’ + f x ' f X ' 
JO Jx • 


T* dx’ 



( 11 ) 
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By reversing the order of integration, the double integral in equation (11) can be trans- 
formed into equivalent single integrals. This is accomplished by first writing the double 
integral in terms of dummy variables of integration 77 and /3 as follows : 



T’frSXOdx' 


dx’ s f x ' f X ' T’(/3,X')d/3 d7? 
JO J-q 


where X' is treated as a constant. Reversing the order of integration leads to 



T’(j3,X’)d/3 d?7 = f x ' ft T’03,X’)di7 d/3 + f X ' f x ' T '{P,X')dr] d/3 

JO JO Jx' JO 

= f x ' pT '(P,X')d(3 + f X ' x'T'(p,X')d0 

Jo Jx' 


so that 


fx * /* x ’ T^x'XOdx’ dx’ = f x ' x , T , (jc , ,X , )<te f + x* f X ' T’(x\X')dx’ (12) 

JO Jx' JO Jx' 


Equation (11) then takes the form 


T’(x\X’) = 1 - - 1 — + 1 


1 + R \dT» 


<L(L f X ' t* dx’ + x’ f X ' T’ dx' + f x 1 x'T' dx') (13) 

1 + R0 t’\R./O Jx' Jo j 


For convenience, let 


I(x’,X’) = — f X ' T’ dx’ + x’ f X ' T' dx’ + f x ' x'T' dx' 

rJo Jx' Jo 


and note that 


01 01 dX' 


0t' 0X' dr' 


As a result, equation (13) takes the form 
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RS 01 dX' 


( 14 ) 


= + 

1 + R \dT’ / 1 + R 3X* dr’ 


When equation (14) is evaluated at the solid-liquid interface x* = X’, there results 

1 = 1 + RX’ /dX’ + \ RS dG(X’) dX* ^ 

1 + R W’ / 1 + “ 


RS dG(X*) dX ' 
R dX’ dr’ 


where 


G(X') ■ I(x* = X’,X’) = - r X ’ T» dx’ + r X ’ x'T’ dx’ 

rJo Jo 


RJO JO 

Solving equation (15) for dX’/dr' gives 

dX’ _ R(1 - X’) 

dT ’ 1 + RX’ + RS dG 


dX’ 


(16) 


Separating the variables and integrating yield 



1 + RX» 
R(1 - X’) 


dX’ + S 


r 

j 0 


1 

1 - X’ 


dG 

dX’ 


dX’ 


(17) 


When the first integral on the right is evaluated and second integral is integrated by 
parts, the final expression for t* = t'(X’) is 


T’ = 


1 + R 


ln(l - X’) 


+ S 


G(X’) 


1 - X’ 



G ( X< ) dX* 

(1 - x *) 2 


(18) 


This is the expression that will be used to compute the frozen layer thickness as a func- 
tion of time. Since G(X’) contains the temperature distribution, this distribution is ob- 
tained from equation (14) and equation (16) is used to eliminate dX'/dr’. The final rela- 
tion for the temperature distribution in the frozen layer becomes 
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(19) 


T’ - R(l - *') | 

1 + R 


1 + Rx’ + RS 


1 + R 


R(1 - X’) 


1 + RX' + RS 


dG 

dX’ 


Equation (19) is an integral equation for the temperature T'(x’,X’). To solve equa- 
tion (18), equation (19) must first be solved so that the integrals that comprise G in 
equation (18) can be evaluated. Two methods for carrying this out are examined herein. 

The first method, and ultimately the more accurate of the two, is an iterative tech- 
nique whereby the integrals in I and G are evaluated by using temperature distributions 
that are obtained from successively better approximations of equation (19). It appears 
that this procedure will ultimately lead to an exact solution. In the second method, equa- 
tion (19) is not used. Instead, the temperature distribution to be used in G(X') of equa- 
tion (18) is approximated by a second-order polynomial with coefficients that are evalua- 
ted using the physical conditions at the boundaries of the frozen layer. Since this approxi- 
mate temperature profile does not arise from the energy equation, the second method 
cannot yield an exact solution. The two methods are outlined in the next sections. 


Solution by Analytical Iterations (AI Method) 


An iterative procedure is utilized by Adams (ref. 10) for some solidification prob- 
lems; before the present method is discussed, a few comments are in order. When 
Adams' procedure was applied to the problem under study herein, the first iteration 
(second approximation) resulted in a cumbersome quadratic equation for dX'/dr’. After 
the quadratic equation was solved, a numerical integration was required to determine X’ 
as a function of t'. A second iteration would have been extremely difficult and was not 
attempted. Hence, an alternate approach was devised in which higher order iterations 
could be analytically obtained. These will now be presented. 

A first-order approximation to the temperature distribution and growth times is con- 
veniently found by neglecting the effect of heat capacity within the frozen layer. When the 
c becomes zero, the parameter S = 0, and equations (18) and (19) reduce to 

Jr 

Tr = -X’ - i-t-5 ln(l - X’) (20) 

1 R 


T , = R(X’ - x*) = RX’ (1 _ ^ 
1 1 + RX' 1 + RX’ 


( 21 ) 
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A second approximation for T” and r' = t'(X') was obtained by substituting equa- 
tion (21) for TJ in equations (18) and (19) while retaining the S terms. After the inte- 
grations in I(x',X’) and G(X’) are performed, the results can be simplified to 


T h = T i - s< 


RX’ L, + R + 2 1 + 3 + 3R + R 2 ln ( X _ x ,^ + ln(RX* + 

3(RX’ + 1)L R(R + 1)J 3( i + R )2 " 3R(1 + F 


i) 

3R(1 + R) 3 , 


( 22 ) 


and 


where 


r l + RX’| + RS 


01 


T , _ R(1 - X'$) 


0X* 


L n 


1 + R 


1 + R 


R(1 - X’) 


1 + RX’ + RS 


dG 

dX’ 


(23) 


dG(X') 


dX’ 


R 


X' | X' 2 t X' 3 


I (1 + RX ') 2 ^ 2 R 


(24) 


91 


R 2 

i| 


9X’ 

I 

(1 + RX’) 2 

r' 

Vr 


[ X» 2 g + X’ 3 A _ [ X ,2 g 2 + X’ 3 | 3 ' 


R 


2R 


(25) 


When the procedure is repeated, that is, when equation (23) is substituted into equations 
(18) and (19), the third approximation results, which is the final approximation evaluated 
herein; that is, 


T in = T i + s 


G n (x») r x ' 

l - x t I 

Jo 


G n (X’) 

(i - xr 


dX’ 


(26) 


1 + RX’| + RS 


T » _ R(1 - X*|) 

™ — mr - 1 


91 

ax’ 


1 + R 


R(1 - X’) 


1 + RX’ + RS 


dG 

dX’ 


(27) 


'n 


where 
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I 


X' 1 R - 1 x» 2 - — 
2 3^ 

1 +R 



X’ _ x ,2 _ RX 

R 3 

1 + R 


X' 2 , 4X' 3 ^ 2X' 4 ^ 2X’ 5 

-f + + 


(1 + R)(1 + RX’) 2 \R 3 3R 2 3R 15 


(x , J = (1 + RX’g)(l - XQ + fw I 

’L 1+R Ut’I, 


(1 + RX’)(1 + RX’g) 
R(1 + R) 


(1 + R)(l + RX’) 2 L\ R 


1 + RX'M / 2X' 5X' 5X' 

R / \ „2 2R 6 


+ /l + RX’Vx’ 2 g 2 , X’ 3 g 3 

\ R /\ 2R 6 


2R 4 S 


(1 + R)(l + RX’) 3 L\ R 


1 + RX'A / _ 5X’ 3 _ 5X’ 4 


,2 6R 24 


+ /x’ | X’ 2 \ /x ,2 | 2 | X’ 3 g 3 \ _ X’ 4 g 4 _ X’ 5 g 5 
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No attempt was made to analytically evaluate the integral in equation (26), since there did 
not seem to be any advantage to doing so because of its algebraic complexity. Since the 
entire equation was to be evaluated on a digital computer, the integration was performed 
numerically. 
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It is obvious from the complexity of the temperature distribution equations (eqs. (27) 
to (31)) that a fourth analytically determined approximation would be almost impossible to 
carry out. As will be shown, a fourth approximate solution is unnecessary because the 
AI method converged rapidly with the first three approximate solutions. 


Approximate Solution by the Prescribed Temperature Method (PT Method) 

The basic expression (eq. (18)) that predicts the dimensionless growth times r' for 
a layer of thickness X' contains integrals in the function G of the instantaneous tem- 
perature distributions in the frozen layer. Because the temperature distributions appear 
as integrands, the growth time r* is not very sensitive to the exact shape of the distri- 
bution. Therefore, using a reasonable approximation for the temperature profile should 
result in fairly accurate predictions for t’ = t’(X'). This idea was first applied to prob- 
lems of freezing and melting by Goodman (ref. 9). 

For the problem being treated in this report, the temperature profile within the 
frozen layer is now approximated by a second-order polynomial of the form 

T’ = A(X’ - x’) + B(X' - x’) 2 (32) 

Equation (32) already satisfies the boundary condition that the temperature should be equal 
to Tf at the liquid-solid interface. Two additional boundary conditions are used to eval- 
uate A and B. The boundary condition at the wall is that the following heat balance be 
satisfied: 
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'x=0 
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which has the dimensionless form 

= r[t*(x' = 0) - l] 

x'=0 

Using equation (32) in equation (33) provides the condition 

-A - 2X’B = R(AX T + BX' 2 - 1) 
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(34) 
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The second boundary condition is at the liquid-solid interface and is given by equation (5) . 
Equation (5) is modified in order to express it entirely in terms of T, without dX/dr, 
by using the following relation that can be written for any location in the frozen layer: 


dT = — 
3x 


dx + H 
3r 


dr 


At the instant the interface moves through position x, x = X, dx/dr becomes dX/dT, 
and the derivative dT/dr = 0 since the interface is always at T^. Then, 
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This is used to eliminate dX/d t from equation (5), so that 
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If the relation 3T/3 t = (k/pc )(3 T/3x ) from the conduction equation is used, this 
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boundary condition can be rearranged into the form 
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In dimensionless form, it becomes 
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X’ 


The temperature profile (eq. (32)) is substituted into equation (36), and the result is 
rearranged into 
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B = — (a.2 — AR_ 
2 \ 1+R 


( 37 ) 


which gives B when A is known. Equations (34) and (37) are solved simultaneously to 
yield A; that is, 
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(38) 


The coefficient A is the root of a quadratic equation, and the positive sign in front of the 
square root was chosen so that A would approach the proper limit as steady state is 
achieved. The limit is found by realizing that the steady-state temperature profile in the 
frozen layer must be linear; hence, from equation (32), B(X' — 1) = 0. Then, from equa- 
tion (37), A(X’ — 1) = R/(l + R), which is also the limit of equation (38). 

Equations (37) and (38) give A and B as explicit functions of X' and, hence, the 
temperature distribution (eq. (32)) is given as a function of X' and x\ The temperature 
distribution is then used to evaluate the integrals in equation (18). 

The quantity G(X') is 


G(X') = x .2 + B x ,3\ AX^ + BX^ 

R\2 3/6 12 


(39) 


while the quantity 
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(40) 


has to be integrated numerically because the A and B contained in G are very com- 
plicated functions of X*. Equation (18) is then evaluated to give the time required to 
form a thickness of frozen layer t* = t’(X') for any values of the parameters R and S. 

With A and B known as a function of X' and, hence, as a function of r', the tem- 
perature profile in the frozen layer can be evaluated at any time during the layer growth 
from equation (32), which can also be placed in the form 


T* = AX’(1 - |) + BX' 2 (1 - |) 2 


(41) 
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Determination of Frozen Layer Growth by Integrating Heat Balance 
at Liquid-Solid Interface (IF Method) 


This alternate approach suggested by Goodman (ref. 9) begins with the interface con- 
dition given by equation (35): 
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This has the dimensionless form 
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The approximate temperature profile (eq. (32)) is used to evaluate the derivatives at the 
interface to give 


dX’ _ 1 + R 2B 
dr’ RS A 

The quantity B is eliminated by using equation (37), and the variables are then separated. 
The result is integrated to yield 


I A - 1 

Jo R 

This was integrated numerically for various R and S using A as a function of X' 
from equation (38). As in the previous method, the temperature distribution is found 
from equation (41) with B obtained from equation (37). 


19 



Comparison of Solutions 


Figure 2 shows the growth of the solidified layer for a few values of the parameters 
R and S as predicted by four different equations presented in the analysis. Included are 
results for the case in which the subcooling in the frozen layer is neglected (S = 0) as 
computed from equation (20). First, the lower set of curves, where R = 0.01, is con- 
sidered. In the limit, when heat capacity is neglected (S = 0), all the solutions obtained 
from equation (18) are in agreement by definition with that given by equation (20). When 
S is increased to a maximum practical value of 3, all the solution methods are in agree- 
ment and fall away only slightly from the S = 0 curve. 

The upper set of curves in figure 2 is for R = 100, and results are shown for the 
complete range of S values. The curves for S = 3 represent the extreme case of a 
large heat-capacity effect and this, coupled with the large R value, provides the maxi- 
mum deviation between the calculation methods. One of the most significant findings is 
the very close agreement between the second and third analytical iterations (curves AIjj 
and AIjjj) . This agreement indicates that the third iteration has converged very close to 
the exact solution, a fact that will be further demonstrated by the temperature distribu- 
tions in figure 3. Hence, for practical purposes, the third iteration can be regarded as 
the correct solution. The method using a prescribed temperature profile (PT curves) is 
in very good agreement with the analytical iterations and, hence, the prescribed temper- 
ature method can also be used with good accuracy for the present problem. The interface 
method (IF curves) deviates somewhat from the other methods but, at a given r% it still 
provides X' values within 5 percent of the correct solution. The reason for this devia- 
tion stems from the fact that the growth times are related to the first and second partial 
derivatives of the temperature at the interface (eq. (42)). For the IF method, these de- 
rivatives are obtained from the approximate temperature profile (eq. (32)) rather than 
from a solution of the energy equation as in the AI method. 


Temperature Distribution Within Freezing Layer 

A few temperature distributions within the freezing layer are shown in figure 3. The 
dimensionless temperature is shown as a function of position within the layer for four dif- 
ferent dimensionless layer thicknesses. An interesting characteristic is that all the tem- 
perature profiles are fairly close to being linear. The heat convected by the warm liquid 
to the interface, combined with the latent heat of fusion that also arises at the interface, 
is conducted through the frozen layer and tends to establish a linear temperature profile. 
The heat capacity of the layer, however, tends to produce a curved profile. The fact that 
most of the temperature profiles are nearly linear suggests that the subcooling energy is 
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small compared with the convective and latent heats at the frozen layer interface. When 
heat capacity is neglected (S = 0), the temperature distribution is linear. 

Figure 3(a) shows the temperature distributions as compared by various methods for 
large values of both R and S. The rapid convergence of the successive analytical iter- 
ation solutions (AI) is demonstrated. The difference between the temperature distribu- 
tions for the second and third iterations is very small. If a fourth approximation were 
evaluated, it appears that it would be indistinguishable from the third approximation and 
would have an insignificant effect on the frozen layer growth curves in figure 2. The pro- 
files evaluated from the PT and the IF methods are not too far removed from the AI pro- 
files. Figures 3(b) and (c) show how the agreement of all the methods is improved as 
either S or R is made smaller. 


Solidification Times as Function of R and S Parameters 

While all the solutions agree well with the third AI, which is taken to be the correct 
solution, they vary considerably in their form and difficulty of evaluation. The simplest 
solution is equation (20) for the S = 0 condition, but its range of application is restricted 
to small R values, as indicated by figure 2. The second AI (eq. (22)) is quite accurate 
for all ranges of R and S and has a form that can be easily and quickly evaluated even 
on a desk calculator if necessary. The third AI, PT, and IF solutions all require a digi- 
tal computer for their evaluation. For these reasons, the second AI solution is the most 
convenient analytical form given here for use in engineering applications. The graphical 
results presented herein were evaluated from the third AI solution (eq. (26)). 

In figure 4 is presented the dimensionless time r’ as a function of R with S as a 
parameter. Each group of curves gives the r’ required to form a given frozen layer 
thickness X'. This set of curves is used later for comparison with the experimental re- 
sults. This can easily be done when there is specified a solidifying liquid and wall ma- 
terial, both having known properties, and the values of a, T^ , T c , h^ , and h c> With this 
information, the quantities R and S can be calculated. From figure 4, the r* values 
can be found by simple interpolation for each X’ shown. Then the relation between di- 
mensional r and X can be calculated from the definitions of r’ and X’ . 

Additional discussion on the general interpretation of the analytical results is given 
in reference 22. 


EXPERIMENTAL INVESTIGATION 
Apparatus 

Flo w circuits . - The test system shown in figure 5(a) consisted of two flow loops, 
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one for water and the other for chilled alcohol. Ordinary tap water was recirculated 
from a 50-gallon insulated storage tank. The thermal capacity of the large water volume 
maintained its temperature constant (hiring a test. The water passed through a flowmeter 
and then a cotton filter to remove small particles that might clog the flow damping screens 
in the test section. 

The alcohol was recirculated through insulated piping from a 50-gallon insulated tank. 
Within this tank was a smaller stainless-steel tank that contained a mixture of dry ice and 
alcohol. By adjusting the liquid level in the inner tank and thereby regulating the extent 
of the cold surface area, the temperature of the alcohol in the outer tank could be adjusted 
to the desired value for a test. 

Test section . - The test section consisted primarily of two rectangular channels that 
were separated by a dividing wall. The important features are shown in figure 5(b). An 
Inconel plate, 4 inches long, 1^ inches wide, and 3/16 inch thick, was embedded in a sheet 
of laminated thermoplastic of the same thickness. This sheet was clamped between two 
channels thereby forming the dividing wall between them (see section BB, fig. 5(b)) . The 
test plate orientation was horizontal. The flow channel walls were made principally from 
1 -inch-thick clear plastic to provide insulation and to facilitate illumination of the ice 
layer that was formed on the Inconel plate . Both sides of the water flow channel adjacent 
to the Inconel plate contained 1^ -inch-thick optical glass windows to minimize distortion 
in the view of the ice layer. A piece of porous filter material and a series of damping 
screens were provided in the water flow channel to produce a uniform velocity profile. 

In early tests, it was found that if care was not taken to have the v^ter velocity uniform, 
the ice layer developed grooves and ridges parallel to the flow direction so that the ice 
thickness was not uniform across the 1^ -inch width of the test plate. With damping 
screens installed, the ice layer was very flat in the transverse direction. The alcohol 
channel, which was not as high as the water channel, contained a piece of filter material 
followed by one damping screen. The uniformity of the alcohol velocity was not as criti- 
cal as for the water because the alcohol was not in direct contact with the ice. 

Instrumentation . - The most important measurements required in the experiment 
were (1) the transient and steady-state temperatures on both surfaces of the Inconel test 
plate, (2) the bulk temperatures of the water and alcohol flowing over the test plate, and 
(3) the ice layer thickness during its transient growth period and at steady state. The 
methods used for each of these measurements are now described in detail. 

Test plate surface temperatures: The surface temperatures of the Inconel plate 
were measured by nine thermocouples on each surface (fig. 5(b)). The pair of 0.003- 
inch-diameter enameled copper -constantan wires forming each thermocouple was cemen- 
ted side by side inside a surface groove 0.007 inch wide, 0.005 inch deep, and 0.75 inch 
long. A thermocouple junction was formed by soldering the wire ends to the plate at the 
end of the groove. Then the surface in the vicinity of the thermocouple junction and along 
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the lead groove was smoothed by removing the excess solder and cement with a very fine 
abrasive. For convenience in handling, the thin lead wires were spliced (using a sol- 
dered joint) to thicker 0.020-inch-diameter wires of the same material. Originally, the 
splices were made in the region between the edge of the Inconel test plate and the outside 
surface of the viewing windows. Preliminary tests revealed that the temperature gradi- 
ent along the length of the splices was causing errors in the surface temperature data be- 
cause of extraneous voltages generated in the soldered splices. For the final test re- 
sults, the splices were made in an isothermal region at room temperature several inches 
away from the outside surface of the test section. 

The method of installing the thermocouples on the plate and the small wire diameter 
were chosen primarily to provide accurate surface temperature measurements during 
both the transient period and at steady state. The small size of the wire allowed the ther- 
mocouple to respond almost instantaneously to surface temperature changes, and mini- 
mized distortions of the local temperature distribution. The length-diameter ratio of the 
wires was made large enough to prevent errors arising from heat leakage from the junc- 
tion through the leads. 

Water and alcohol bulk temperatures: The bulk temperatures of the water and alco- 
hol were measured with thermocouple probes at the inlets of their respective flow chan- 
nels in the test section. The alcohol temperature was also checked with a thermocouple 
at the channel outlet. In the storage tanks, thermometers provided the temperature of 
each fluid as an approximate check on the values in the test section. The probes and 
thermometers were calibrated in a water-crushed ice mixture and at lower temperatures 
down to -60° F by comparison with a precision thermometer in a dewar of cold alcohol. 
All thermocouples could be read with either a precision hand potentiometer or a strip 
chart electronic recorder. 

Ice-layer thickness measurements: Three basic methods were employed to measure 
the ice-layer thickness during the transient growth period and at steady state. These 
will be referred to as the (1) photographic (2) visual, and (3) temperature probe methods. 
These methods will be discussed in some detail for the benefit of those who may want to 
make similar measurements. 

The photographic method that was used for some preliminary results employed a 
70-millimeter camera equipped with a timer, automatic film advance, and strobe lights 
(fig. 6). With this equipment, a series of still photographs could be taken at predeter- 
mined time intervals throughout the growth period. The camera viewed the ice layer 
edgewise, as shown in figure 6. The photographs were enlarged to 8- by 10-inch size, 
and the ice thickness measured from them. This method proved to be of limited accuracy 
because of optical distortion effects, optical illusions, difficulties in obtaining exact cam- 
era alinement with the plane of the ice interface, and limited image magnification. It is 
worth mentioning that a 16-millimeter motion picture camera was tried but was unsatis- 
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factory because the pictures were too small to provide the necessary clarity. Generally, 
it was found easier and more accurate to obtain results by the visual method and with the 
temperature probe described next. 

The visual method of measuring the local ice -layer thickness utilized a cathetometer 
and two film sheets with a grid printed on each (fig. 6) . The grids consisted of a group 
of carefully spaced parallel horizontal reference lines with other lines at a 45° and 90° 
angle to those reference lines. The films were photographic reductions of a grid that was 
carefully ruled on a large piece of white poster board. The horizontal reference lines 
were 0. 025 inch apart and were coded for easy identification of their position relative to 
a zero reference line. Short horizontal dashes spaced vertically 0.005 inch apart were 
made between the horizontal reference lines at 1/2 -inch intervals along them. These 
grids were mounted on the outside of the optical glass windows and were positioned so 
that the zero reference line on each grid was in the plane of the top surface of the Inconel 
test plate. Hence, when the grids were viewed with a cathetometer telescope, two cor- 
responding reference lines, one from each grid, were in the same plane at the correct 
height above the Inconel plate. To measure the ice thickness at any position along the 
layer length, the ice -water interface was sighted with the telescope and the interface po- 
sition was read from the film grids. 

In actual practice, two factors affect the accuracy with which the visual (and also the 
photographic) method can be used to measure ice-layer thicknesses. These were espec- 
ially significant for measurements made during the transient growth period and for layers 
less than about 0. 050 inch thick. The factors are the transparency of the ice layer and 
the optical distortion that resulted from the small density gradient in the thermal boundary 
layer above the ice-water interface. The 45° and 90° lines on the film grids were for the 
purpose of assessing the extent of any optical distortion that might be present. 

The transparency of the layer made it difficult at times to detect the exact location 
of the ice -water interface. The interface definition was aided by judicially arranging the 
lighting and using a white background. It was found that the light positions had to be 
shifted for various tests. When the test conditions were changed, it was difficult to pre- 
dict exactly how to adjust the photographic lighting to assure sharp interface definitions 
throughout the transient. This was a drawback in the photographic method since, after 
the test had been completed and the films developed, some of the images were sometimes 
lacking in clarity, and a complete set of transient data could not be measured. For this 
reason, the use of the cathetometer and film grids offered a distinct advantage over the 
use of photographs. 

The optical distortion resulting from bending of light by the nonisothermal water dis- 
placed the image of the interface from its true position. To account for optical distortion 
and to have the means for checking and calibrating the visual thickness measurements, 
several devices were employed. A few reference lines at known heights above the Inconel 
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plate were scribed on the inside of the optical windows. Some metal pins of known heights 
were mounted on the plastic plate (that contained the Inconel inserts) downstream of the 
ice layer. These provided reference lengths inside the test section adjacent to the ice 
layer. A piece of stainless-steel rule with 1/64-inch divisions was also mounted verti- 
cally on the plastic plate to the rear of the Inconel test insert. For measurements of 
steady-state ice thicknesses, a micrometer head was mounted on top of the water chan- 
nel, and a vertical traversing rod was extended down into the water stream. The rod po- 
sition was above one of the thermocouple locations in the center region of the Inconel 
plate . After an ice layer had formed, the rod that had been initially withdrawn from the 
water stream could be moved down to touch the ice surface and the thickness could be 
measured with the micrometer. The rod had to be moved very quickly since it would dis- 
tort the water flow pattern near the ice surface and would also conduct heat toward the 
interface; both these actions would tend to melt the ice away before the probe could touch 
it and the thickness reading could be obtained. 

The optical distortion caused the grid lines that were viewed through the water to 
bend significantly in the vicinity of the ice -layer surface. To calibrate the system, a 
steady ice thickness was formed, and the cathetometer was sighted across the interface 
at the position along the layer where the micrometer probe was located. Thickness read- 
ings were then taken from the two grids and with the micrometer (a thermocouple probe 
described below was also used as a further check) . It was generally found that the grid 
viewed through the water gave the thickness accurately while the grid on the front window 
read about 0. 020 inch too large. 

The third method used to measure the ice-layer thickness employed a thermocouple 
probe (fig. 6). The probe consisted of a butt -welded junction made from 0. 003-inch 
Chromel constantan wires stretched between two needles mounted on a 0. 100-inch- 
diameter shaft. The shaft extended through the top of the test section and was mounted in 
a micrometer head, so that the junction could be accurately positioned at any specified 
height above the test plate . By touching the ice layer with the junction, the layer thick- 
ness at that location could be determined. The section Test Procedure describes the 
measurement of transient ice thicknesses by use of the thermocouple probe. 

The temperature probe was also used to check the readings of one of the Inconel plate 
surface thermocouples. This was done by lowering the probe junction to the surface; by 
design, the probe location was very close to one of the surface thermocouples. 


Test Procedure 

O btaining transient data . - All the transients were initiated by a sudden increase of 
the alcohol flow. For most of the tests, a very small alcohol flow was initially present 
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to provide a small, very thin patch of ice on the plate. For a few tests, there was no ini- 
tial alcohol flow and no ice patch; on starting the alcohol flow, the plate subcooled below 
the freezing point and ice suddenly flashed on. 

The transient data was taken only in the central region of the plate length where the 
heat transfer was essentially one dimensional in conformity with the conditions of the 
analysis. The heat transfer through the ice layer and Inconel plate was two dimensional 
near both plate ends. 

For a measurement of the instantaneous thickness, the cathetometer and grid films 
were first alined. The cathetometer was sighted at one of the thermocouple stations on 
the plate. The full alcohol flow was then turned on, and the moving ice interface was 
viewed with the telescope. The time when the ice interface moved past each grid line was 
noted, and the thickness was marked at that time on the temperature strip chart that re- 
corded the upper and lower plate surface temperatures. Thus, the strip chart provided a 
record of the thickness and plate surface temperature variations with time. 

An alternate means for measuring the transient thicknesses was by use of the thermo- 
couple probe shown in figure 6. This method served as a check on the data obtained with 
the cathetometer. The use of the probe was possible for two reasons: (1) the thermo- 
couple and its needle -like supports were sufficiently small so that no local distortions of 
the interface were observed when the thermocouple was placed in contact with it and 
(2) the interface temperature measured with the probe was always found to be between 
31.7° and 32° F for transient and steady-state conditions. To measure the instantaneous 
thickness during the growth period, the probe junction was positioned 0. 005 inch above 
the plate, and the probe temperature was monitored with a fast response recorder. As 
the ice layer grew, the probe temperature decreased and reached 32° F when the probe 
was contacted by the ice. The probe was then quickly raised 0. 005 inch before it could 
freeze into the advancing layer. A step jump was produced on the recorder as the probe 
moved upward into the warm water stream. This procedure was repeated giving a series 
of steps on the recording chart at times corresponding to every 0. 005 inch of the growth. 
The temperature probe method proved to be the best method for measuring the instan- 
taneous thickness during the transient growth period. The measured thicknesses were 
accurate to within 0. 001 inch; they were easy to obtain, and agreed well with those mea- 
sured with the grids and telescope. Because the probe was divised relatively late in the 
duration of the experimental project, most of the data reported here were obtained with 
the cathetometer. 

Calibration for heat l eakage . - Although the sides and windows of the test section 
channels were about 1 inch thick, some heat leakage was expected from the room into the 
cold test plate and alcohol. Of even greater importance was the leakage from the warm 
water through the walls of the test section and into the side edges of the Inconel test plate. 
As a result, some of the heat being removed by the coolant originated from heat leakage 
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and not from the ice formation process. To determine a correction for heat leakage, the 
following procedure was utilized. A piece of polyurethane insulation was formed into the 
shape of an ice layer and glued to the surface of the Inconel test plate on the water side. 

In this way, the flow pattern in the water stream was comparable to that during a test with 
an ice layer, and the direct heat flow through the test plate was small. With this arrange- 
ment, the plate was found to be from 0° to 6° F above the alcohol temperature depending 
on the temperature levels of the alcohol and water; without heat leakage, the plate would 
be at the alcohol temperature. The heat leakage temperature difference between plate 
and alcohol was measured for various flow rates and various alcohol and water tempera- 
tures. As shown in appendix B, it is necessary to add this temperature difference to the 
measured alcohol temperature to obtain the alcohol temperature that is effective in the 
ice growing process. 


Processing of Data 

To compute the expected ice growth from theory for comparison with experiment, the 
parameters R and S must be known along with the quantities that relate dimensionless 
to actual thicknesses and times. Four important quantities are the thermal conductivities 
of the Inconel plate and the ice, and the heat-transfer coefficients on the water and alcohol 
sides. 

Thermal conductivities. - Data for the Inconel plate thermal conductivity were ob- 
tained from references 23 and 24. For the range of the present tests, the plate temper- 
atures varied from 32° F at the initiation of a transient to about -20° F at the end of a 
transient. For this temperature interval, the Inconel thermal conductivity k w was taken 
as 7. 8 Btu per hour per foot per °F. 

The thermal conductivity of the ice k was obtained from Ratcliffe (ref. 25) and was 
selected as 1. 36 Btu per hour per foot per °F. This was chosen as an average value for 
the temperature range 32° to -10° F ordinarily encountered in the present test series. 
Some comments regarding the ice conductivity are given later in connection with the test 
results. 

Heat-transfer coefficients . - To compare the measured ice growth with the theo- 
retically predicted values, the local convective heat-transfer coefficients on both the 
water and alcohol sides must be known for each flow rate. The coefficients must be known 
with good accuracy, since the predicted ice growth can be quite sensitive to them. The 
sensitivity arises from the form of equation (1), which shows that the steady-state thick- 
ness depends on the difference of two terms. As this difference becomes small compared 

with the magnitudes of the terms, the X value becomes increasingly sensitive to the h ; 

s c 

and h values, 
c 

A first attempt to determine h^ and h c values was by using laminar boundary layer 
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theory for flow over a flat plate of uniform temperature. This always produced theo- 
retical values about 25 percent below those measured experimentally by methods that are 
described next. The discrepancy could have resulted from a number of experimental 
factors such as imperfections in the uniformity of the experimental velocity profile, free 
convection, disturbances from the side walls of the flow channels, variable fluid proper- 
ties, etc. An additional complication was that the cold alcohol absorbed air, and some 
bubbles would be released when the alcohol contacted the warm plate. This would in- 
crease the heat -transfer coefficient above the laminar value. For these reasons, it was 
thought best to experimentally calibrate the apparatus to determine the convective coef- 
ficients. 

Two methods were used for experimentally obtaining the h^ and h c values. One 
approach was to set the flow rates for both water and alcohol at values to be used in the 
transient ice formation tests but not have the alcohol cold enough to form ice. Then, the 
temperature drop through the Inconel plate was measured at the locations along the plate 
length where the ice growth would later be observed. The known plate conductivity was 
used to compute the heat flow through the plate, assuming one -dimensional heat flow, and 
the local heat -transfer coefficients on both sides of the plate were found by using the mea- 
sured temperatures of the liquids. A shortcoming is that the h^ values determined by 
this method are for a nonuniform surface temperature condition, whereas the boundary 
condition with an ice layer is one of the uniform temperature. Consequently, the h^ 
values from a bare plate determination are only approximately applicable to conditions 
with an ice layer. 

A second approach was to grow steady-state ice layers of various thicknesses by 
holding the flow rates constant and allowing the water or alcohol temperatures to be var- 
ied from one test to another. Both the ice layer and Inconel plate thicknesses could then 
be used as heatmeters to obtain local heat fluxes and the corresponding local convective 
coefficients. For a given water and alcohol flow rate, usually three different ice thick- 
nesses were formed by, for example, holding the alcohol temperature constant and using 
three different water temperatures. The corresponding heat -transfer coefficients had 
about a 10-percent spread in values, and an average was taken to provide the coefficients 
characteristic of the fixed flow rates. 

For the same flow rate, the deviations between the coefficients determined by the 
second method and those using a bare plate were a maximum of 12 percent. There was 
no criterion by which either of the two methods could be shown to be the best. Hence, 
for a given position on the plate, the local values determined from the two methods were 
averaged and used in the theoretical predictions. It should be noted that the effective 
coolant temperature as described in appendix B was the coolant temperature used for all 
processing of data for comparing experimental results with theory. 

Use of data in theore tical predictions . - The theoretical predictions of X, T^, and 
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Tg involved the following important assumptions: (1) the h^ and h c remain constant 
throughout the transient period, (2) the ice-liquid interface is always at 32° F, and 
(3) there is no contact resistance between the ice and the Inconel plate. Some discussion 
on the validity of these assumptions is given later. With a given set of test conditions h^ , 
h c , , and T c , the quantities X g , X(r), T^r), and T 2 (r) can be computed from the 
analysis. The theoretical steady-state ice thickness is computed from equation (1). The 
R and S values can then be found. Figure 4 is used with these R and S values to find 
the dimensionless time r* corresponding to various X’ values. The theoretical X(t) 
variation is then found by using the definitions of the dimensionless variables to convert 
back to dimensional variables. 

The calculation of the theoretical variations of Tj(t) and T 2 (t) was simplified be- 
cause, for the range of experimental conditions encountered in the present series of 
tests, the subcooling parameter S was always less than 0. 25. For this small subcool- 
ing and because R was always of order unity, figure 3 shows that the temperature T 2 
at the interface between the ice and plate was always very close to the values for the 
limiting solution when S = 0. For S = 0 and x’ = 0, the temperature equation (eq. (19)) 
simplifies to 


T' = T2 " Tf = R _ / 1 \ R(1 - X’) 

T c - T f 1 + R \1 + R/ 1 + R X’ 

This yields the expression for T 2 

To = T f + (T - T f ) - P*' (44) 

1 1 c 1 1 + RX' 

that was used for comparison with the experimental temperature trace at the ice-plate 
interface . 

To calculate the theoretical Inconel-alcohol interface temperature T^(t), the heat 
conduction through the Inconel plate must be examined. One basic assumption made in 
the Analysis section was that the plate heat capacity was negligible. The validity of this 
assumption was examined for typical data by using equation (C 5) developed in appendix C . 
The term 


P w C P>w _ 9_ ( 2T + T ) 
6 0r 1 £ 


in equation (C5) was evaluated by using the measured Tj(t) and T 2 (t) shown in figures 
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8(a) and (b) and one of these results (for fig. 8(b)) is presented in figure 10. The maxi- 
mum heat rate extracted to subcool the plate was less than 10 percent of the total heat 
flux through the plate. This maximum occurred early during the transient, and the sub- 
cooling contribution then diminished rapidly with time. The maximum correction to the 
calculated T^(t) was less than 1. 5° F. For this reason, the values for Tj(t) shown in 
the figures were calculated from a simplified form of equation (C5) in which the heat 
capacity term was neglected, that is. 


h (T, -TJ=-W(T 9 -T 1 ) 


Therefore, solving for Tj, 


T 


1 ~ 



+ h T 
c 


c 



(45a) 


(45b) 


RESULTS AND DISCUSSION 
Configuration of Steady-State Ice Layer 

Before discussing the comparison of experimental data with theory, it is of interest 
to consider the shape of the steady-state ice layer that formed on the chilled plate. The 
shape is illustrated in figure 6(b) for the case where the entire plate is below the freezing 
temperature. The convective heat-transfer coefficient of the warm water is very high 
near the leading edge of the plate as is characteristic of flat -plate boundary-layer heat 
transfer. This causes the ice layer to have zero thickness very close to the leading edge. 
If the convective transfer from the warm water is high enough to raise a portion of the 
plate above freezing, then the leading edge of the ice layer is downstream from the plate 
leading edge. Since the convective heat transfer decreases with length along the plate, 
the ice thickness increases in the downstream direction. The ice -layer shape ontour 
resembles that of a flow boundary layer (similar contours were obtained by Zerkle'in 
ref. 16). The ice layer extended a short distance past the trailing edge of the plate, as 
shown in figure 5(b). The data in the present experiments were taken in the central re- 
gion of the Inconel plate where the ice growth was very nearly one dimensional. 
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Transient Ice Growth and Surface Temperatures 


Several typical sets of experimental data are now compared with theoretical predic- 
tions. The ice -layer thickness and surface temperatures of the Inconel plate are given as 
a function of time. Zero time is always taken at the moment when ice first appears at the 
location on the plate where the transient measurements are being taken. 

The analysis assumes that, at zero time, the surface of the plate in contact with the 
water is at the freezing temperature and ice is beginning to form. In order to have the 
experimental conditions correspond as closely as possible with the assumptions of the 
analysis, it was necessary to provide a nucleus of ice on the plate by using a very low 
alcohol flow so that, when the full alcohol flow was started, the ice layer would imme- 
diately begin to spread over the plate. If this were not done, the plate would subcool as 
much as 12° F before ice would suddenly flash on the plate. The flashing of ice was 
accompanied by a step rise in plate temperature that was never large enough to raise the 
plate surface to 32° F. With a small patch of ice initially on the plate, the surface would 
still subcool about 1° F while the ice was spreading over the plate, but the initial condi- 
tion duplicated within reason that assumed in the analysis. All the runs given herein, 
with one exception, are for the condition where a small ice nucleus was initially provided. 

The results in figures 7(a) to (c) and 8(a) to (c) are for three test runs with approxi- 
mately the same water and coolant temperatures and coolant heat-transfer coefficient, 
but with three different values of the heat-transfer coefficient in the flowing water. Fig- 
ures 7(a) to (c) give the transient ice thickness, while figures 8(a) to (c) provides the cor- 
responding variations of the plate surface temperatures. The h^ was varied by changing 
the water velocity, which had mean values of 0. 18, 0. 26, and 0. 35 feet per second for 
figures 7(a) and 8(a), 7(b) and 8(b), and 7(c) and 8(c), respectively. 

An inspection of figures 7(a) to (c) indicates a generally favorably agreement of ex- 
perimental ice thickness with theory and, at any time, the theory would provide a suitable 
thickness estimate for engineering purposes. At this time, only the theoretical curve 
given by the solid line is being considered; this line was predicted as described in the 
Analysis. The dotted prediction was computed in part from the experimental data and is 
discussed later. All the experimental data exhibit the same characteristics when com- 
pared with theory (solid curve) . The experimental growth starts a little more slowly than 
predicted and, when the thickness is about 15 percent of the steady-state value, the growth 
becomes more rapid than predicted. Sometimes, as in figure 7(b), the steady thickness 
is achieved much more abruptly than the gradual asymptotic approach that theory indi- 
cates. As expected, figure 7(a) with a low h^ provides a thicker steady-state layer than 
figure 7(c), which corresponds to a higher and consequently has more heat convected 
to the ice -water interface. 

Figures 8(a) to (c) show that the plate surface temperatures in all cases fall below the 
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predicted curves. The difference is generally of the order of 4° F. That the plate tem- 
peratures were colder than predicted from theory is consistent with the observed more 
rapid ice growth rates. Possible reasons for these discrepancies are discussed later. 

Figures 7(d) and 8(d) show some data where the coolant temperature was not as low 
as that used for the previous figures. This provides a smaller steady-state ice thick- 
ness, which results in an R value about one-fourth of that in figures 7(a) and 8(a). The 
comparison with theory is of the same nature as in the previous figures, that is, more 
rapid experimental growth rates than theory with a tendency toward abruptly reaching a 
steady state, and experimental surface temperatures usually about 4° F below theory. 

Figures 7(e) and (f) and 8(e) and (f) show two sets of data that were obtained for prac- 
tically identical values of heat -transfer coefficients and liquid temperatures. The only 
difference is that the test in figures 7 (e) and 8(e) was started without having an initial 
patch of ice on the surface. There is little difference in the transient ice growth, but 
there are some differences in the surface temperatures. This is shown in additional de- 
tail by the small inserts in figures 8(e) and (f). 

The inserts show the measured temperature history for a short time before and after 
the ice forms at the measuring location. Figure 8(f) for the case with an initial ice patch 
shows that, after the full alcohol flow is started, the lower surface temperature T^ de- 
creases monotonically. The upper surface (temperature Tg) subcools to 30. 5° F and 
then remains at that temperature for a few seconds. The temperature Tg then decreases 
as the ice layer is growing in thickness. Figure 8(e), where there is no initial ice patch, 
exhibits a different behavior. The upper surface of the plate subcools to 20° F, and ice 
is then nucleated. The temperature Tg quickly rises several degrees and then begins to 
drop off as the ice layer thickens. The lower surface of the test plate also experiences a 
small temperature rise at the time that the ice is nucleated. For the transient without an 
initial ice patch, the surface temperature throughout most of the transient remains about 
2° to 4° F further below the theory than the temperatures obtained with an initial ice 
patch. 

With regard to the sudden nucleation of ice without an initial ice patch, some inter- 
esting supplementary temperature data are shown in figure 9. The temperature -time 
traces shown in figure 9(a) were made by having the thermocouple probe within the ther- 
mal boundary layer of the water during a time span from shortly before the ice formed 
until after the probe became embedded in the ice layer. The measurements were made 
as follows: The cold alcohol flow was reduced so that no ice existed anywhere on the 
plate, and the probe was positioned at the desired location close to the plate surface. The 
full alcohol flow was suddenly introduced, and the probe output was continuously moni- 
tored on a strip chart recorder . Since only one thermocouple probe was used, the traces 
shown in figure 9(a) for various positions above the plate were taken in successive tests, 
but the temperature magnitudes c?.n be compared because the tests were quite repro- 
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ducible. The time origin t = 0 in figure 9(a) was chosen to be at the local minimum of 
each trace since the abrupt temperature rise always occurred very close to the instant 
that the ice was observed to flash on the surface. Also shown is a cross plot of tempera- 
ture against distance from the wall surface at r = 0 and r = 2 seconds (fig. 9(b)). 

The results reveal the existence of a layer of subcooled water in contact with the wall 
prior to ice nucleation and a subsequent instantaneous temperature increase in the layer 
coincident with the sudden appearance of the ice layer. This temperature increase was 
probably caused by the fact that ice crystals are known to propagate quickly through sub- 
cooled water (ref. 19), and the heat of fusion that is released raises the mixture of water 
and ice crystals to the freezing temperature. A transient then follows in which some of 
the heat flow within the water boundary layer serves to adjust the temperature profile to 
the suddenly increased boundary temperature. 


Discussion of Discrepancy Between Theory and Experiment 

Various factors were investigated to try to determine why a deviation existed between 
theory and experiment and, in particular, why the plate temperatures were consistently 
below the theoretical values. These factors and some interesting supplemental results are 
now discussed. 

Contact resistance betw een ice and plate . - A contact resistance between the ice and 
metal plate could account for part of the temperature discrepancy. Such a resistance was 
noted in the thermal conductivity experiments of Ratcliffe (ref. 25) when metal plates 
were bonded to the ends of thermal conductivity specimens. Several tests with the pres- 
ent apparatus were made to detect the presence of a contact resistance. Before forming 
any ice, the 0. 003 -inch-diameter thermocouple probe was positioned about 0. 002 inch 
above the Inconel surface at a location along the plate close to one of the plate surface 
thermocouples. An ice layer was then formed, and the thermocouple became encased 
within the layer. For steady conditions when the layer had reached its final thickness, 
the probe temperature was within 0.7° of the surface temperature Tg. During the early 
period of ice layer growth, the probe temperature was 2° to 3° F higher than the instan- 
taneous Tg. The fact that the probe was 0. 002 inch above the plate would account for a 
substantial portion of this temperature difference; hence, the effect of contact resistance 
appears small enough to be neglected in the present study. 

Temperature at moving and steady-stat e ice interface . - In the analysis, it was as- 
sumed that the moving ice interface was always at the equilibrium freezing temperature 
(32° F). Conceivably, for a very fast moving interface, the heat removal rate could be 
so large that the crystallization process would not be rapid enough to prevent the interface 
from becoming subcooled. No evidence of appreciable interface subcooling was detected 
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in the present experiments . Tests were conducted by setting the thermocouple probe at 
various distances from the cold surface and then growing an ice layer. The probe tem- 
perature was continuously monitored on a strip chart recorder and was always found to be 
within a few tenths of 1° F of 32° F at the instant that the interface was observed to come 
into contact with the probe. It is interesting to note that, as the ice contacted and grew 
around the probe, the temperature remained constant at about 31.7° F for 1/2 to 2^ sec- 
onds depending on the cooling rate. This indicates that, adjoining the interface, there 
probably exists a thin region in which the crystallization process is occurring. The as- 
sumption in the analysis of an equilibrium temperature at the moving interface appears to 
be satisfactory. 

The temperature of a stationary ice interface in flowing warm water was checked in 
two ways. By using a reduced alcohol flow, a small patch of ice was formed on the test 
plate such that the edge of the patch touched one of the surface thermocouples. For this 
condition, the thermocouple always indicated 32° F within the accuracy of the thermo- 
couple calibration. A second approach was to touch the steady-state ice layer with the 
thermocouple probe. This approach always yielded a reading that was within 32±0. 2° F. 

Transient variation in heat -transfer coefficients . - A possible cause of experimental 
deviations from theory is that the convective heat -transfer coefficients, assumed constant 
in the analysis, are somewhat variable during an experimental transient. As the ice 
layer grows, the temperature of the Inconel test plate decreases and, consequently, the 
average temperature decreases in the alcohol thermal boundary layer. This produces an 
increase in alcohol viscosity and a decrease in coolant heat -transfer coefficient as the 
transient proceeds. 

An additional experimental complication was the absorption of air by the cold alcohol. 
In the early part of the transient when the Inconel plate was still relatively warm, the 
contact of the alcohol with it would cause some bubbles to be released, thereby producing 
turbulence. The bubble release would decrease as the transient proceeded and the plate 
temperature decreased. This would enhance the alcohol coefficient in the early portion 
of the transient. 

On the water side, the ice layer growth produces a reduction in the flow cross sec- 
tion that is a maximum of 20 percent corresponding to the thickest ice layer tested. This 
would increase the water velocity and the convective coefficient during the transient. 

Both a higher alcohol-side coefficient and a lower water-side coefficient during the ear- 
lier portion of the transient would cause a relatively more rapid ice growth during this 
part of the transient. 

Heat fluxes . - A further comparison of measured and predicted results was made by 
examining the heat fluxes through the ice layer and wall and into the coolant. These 
fluxes are exactly equal at steady-state conditions only. During the transient period, the 
heat rates required to subcool the ice and the wall are also present. The heat leaving the 
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ice layer and entering the wall includes the convective heat h^ (T^ - Tj), the heat of fusion 
liberated at the ice-water interface, and the heat for subcooling the ice. Similarly, the 
energy transferred to the coolant h c (T^ - T c ) consists of the heat entering the wall plus 
the energy to subcool the plate. When the subcooling of the ice and plate is small, which 
was true throughout most of the transient period, the instantaneous heat flux through each 
can be computed from the instantaneous temperatures inserted into the steady-state con- 
duction equation. 

Using this type of quasi -steady -state calculation, a typical set of heat fluxes was cal- 
culated from the experimental data given in figures 7 (b) and 8(b) . The fluxes are shown in 
figure 10 along with the theoretical curve. The data in other parts of figures 7 and 8 
yielded similar heat flow curves. For the major portion of the growth period (following 
about the first 80 sec), the instantaneous heat flux through the wall agrees with the flux 
into the coolant as it should. These portions of the curves, however, are displaced to 
earlier times than the theory; this is consistent with a similar behavior for the surface 
temperatures and the more rapid observed ice growth. The explanation for this dis- 
placement is not known. For. times earlier than about 80 seconds, the heat flow through 
the wall is larger than that to the coolant, indicating that possibly a larger h c should 
have been used to compute q during the early portion of the transient. The larger h 
could have been caused by the air bubbles leaving the alcohol adjacent to the plate as men- 
tioned earlier. The bubble release was larger during the early portion of the transient 
when the plate temperature was the highest. 

The largest discrepancy in figure 10 is the computed heat flux through the ice, which 
is larger than the other heat fluxes, especially at early times. This indicates that per- 
haps the ice conductivity value used to compute q^ is too large. A smaller k would 
provide an extra heat flow resistance across the ice layer, which would account for the 
low surface temperatures measured but would not account for the more rapid experimen- 
tal ice growth. There is some foundation for the belief that k would be low during early 
ice growth. Early during the transient when the ice growth was rapid, the layer contained 
some air bubbles, and the ice and its surface appeared to be diffuse. As the layer be- 
came thicker and the growth was slower, the ice was clear and its surface was very 
smooth and transparent. In the ice thermal conductivity tests of Powell (ref. 26), low 
conductivities were detected when ice was initially frozen on a plate. It was found nec- 
essary to apply heat to the plate to melt the ice adjacent to the plate and then refreeze it 
in order to obtain conductivities consistent with those normally quoted in the literature. 
Thus, it is suggested that the ice with the diffuse appearance has a different k than the 
clear ice. 

Ice growth computed from measured surface temperature. - To provide an additional 
examination of the consistency of the experimental data, computations were carried out 
where the ice growth was computed using the measured surface temperature Tg(r) as a 
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known boundary condition. This calculation eliminates the uncertainties involved in the 
Inconel plate conductivity and the alcohol heat-transfer coefficient and determines whether 
the observed ice growth is in accord with a theory using a variable surface temperature 
as the boundary condition. In the present computations, the ice heat capacity was neglec- 
ted because the S values were small for the present experiments. The analytical ex- 
pression used was obtained by evaluating equation (7) at x = X and omitting the heat 
capacity integral on the right side. This gives, after rearrangement, 


dX 

dr 




(46) 


This equation was integrated numerically using the experimental Tg variation. The re- 
sults are shown as a dot-dashed line in figure 7 . The agreement with the data is generally 
as good and sometimes a little better than the prediction using the coolant heat-transfer 
coefficient. The observed ice thickness at any time is always less than what would be ex- 
pected from the recorded surface temperature variation. This indicates, as was pointed 
out previously, that, during part of the transient period, the ice conductivity may be lower 
than the conductivity values reported in the literature . It is evident that there are sev- 
eral factors in the experiment that could, in a combined way, account for the observed 
deviations from theory. To better understand the details of the freezing process, more 
detailed measurements would be desirable, possibly with a thermocouple rake to obtain 
the details of the temperature variations within the ice. The range of experimental op- 
erating conditions should also be extended, for example, to lower coolant temperatures. 


CONCLUSIONS 

Analysis 

The goal of the analysis was to develop a means for predicting the transient growth 
of the frozen layer that forms when a flowing warm liquid is in contact with a cold flat 
plate that is convectively cooled on the opposite side. This goal was achieved by employ- 
ing the following two basic approaches: (1) by integrating the transient conduction equa- 
tion over the entire frozen layer thickness and (2) a heat balance at the frozen layer - 
liquid interface. 

By use of a proper dimensionless layer thickness and time variable, the dimension- 
less frozen layer growth can be expressed as a function of two parameters. One param- 
eter provides a measure of the maximum possible subcooling energy of the frozen layer 
as compared with the latent heat of fusion. The second parameter is the ratio of the heat 
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flow resistance of the steady-state layer to the combined resistance of the wall and con- 
vection coefficient on the coolant side. 

The results of the analysis led to a number of conclusions. Two of the five solutions 
that were derived were the most convenient to use and the most precise. Both were suc- 
cessive iterative approximations that were analytically derived. One of these approxima- 
tions was sufficiently simple that it could easily be evaluated on a desk calculator, while 
the other appeared to converge to the correct solution and was used to prepare the final 
graphical results. These results suggest that the analytical iterative technique developed 
herein is a rapidly converging means for solving nonlinear transient freezing problems. 

There are two additional conclusions that can be drawn from the analytical results. 
The first is that the subcooling of a frozen layer during its formation slows the growth 
rate substantially 'for some conditions. This is particulary true of thick layers. The 
second is that the temperature profiles in the frozen layer were found to be almost linear 
at all times. This suggests that the energy required to subcool the frozen layer is much 
less than the combined convective and latent heat passing through the layer. 


Experimental Study 

In the experimental study, ice was formed on an Inconel plate cooled by flowing 
chilled alcohol. Measurements were made of the instantaneous ice -layer -thickness and 
plate surface temperatures throughout the transient growth period. The most significant 
findings are summarized as follows : 

1. The experimental transient ice thicknesses any any time during the ice growth 
were generally within 15 percent of the theoretically predicted values. However, the 
following departures between theory and experiment were evident: 

a. The experimental ice growth was generally more rapid than predicted, and 
the steady-state frozen thickness was reached more abruptly than expected from the theo- 
retical asymptotic approach to steady state. 

b. The measured temperatures on both surfaces of the Inconel test plate were 
about 4° F below theoretical predictions throughout most of the transient. 

2. Very little contact resistance was found at the interface where the ice was frozen 
onto the Inconel plate. 

3. The surface of the ice for both growing and steady-state ice layers in flowing 
warm water was always found to be, within a fraction of a degree, at the equilibrium 
freezing temperature 32° F. 

4. If there is no ice initially on the metal plate, the plate will usually subcool, as 
much as 12° F in some instances, before ice is nucleated. Ice will then rapidly flash 
over the entire plate. This is accompanied by a sudden increase in surface temperature, 
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but the surface temperature never rises to the freezing point. Following the initiation of 
ice in this manner, throughout most of the transient, the temperatures remained further 
below predicted values than the measurements taken with an initial ice patch. 

5. During the early stages when the ice growth was rapid, the layer contained some 
air bubbles and the ice surface appeared diffuse. As growth became slower, the layer 
was very clear and the surface very smooth and transparent. 

6 . Care must be taken when making visual measurements of ice thicknesses in flow- 
ing warm water because of optical distortion near the ice interface. 

7. A computation of instantaneous heat flow through the ice indicates that the ice con- 
ductivity during the early growth period may have been smaller than later in the transient. 
This is also indicated by the ice thickness variation computed from the measured surface 
temperature. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 12, 1966, 

129-01-09-02-22. 
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APPENDIX A 
SYMBOLS 


A coefficient of linear term in ap- 
proximate temperature profile 

a thickness of the wall on which 

frozen layer is forming 

B coefficient of quadratic term in 
approximate temperature pro- 
file 

Cp specific heat at constant pressure 
of solidified material 

Cp w specific heat of wall on which 
frozen layer is forming 

h convective heat -transfer coef- 

ficient 

k thermal conductivity of solidi- 

fied material 

thermal conductivity of wall on 
which frozen layer is forming 

L latent heat of fusion for solidi- 
fying material 

q heat flow per unit area in the 

positive x direction 

R dimensionless thermal 

resistance parameter, 

(X g /k)/ [(l/h c ) + (a/k w )] 

S dimensionless subcooling param- 

eter, c p (T f - T c )/L 

T temperature 

T* dimensionless temperature, 

(T - T f )/(T C - T f ) 


T c coolant temperature (effective 

value for solidification) 

T c measured (uncorrected for heat 

leakage) coolant temperature 

v fluid velocity 

X thickness of frozen layer 

X’ dimensionless thickness of frozen 

layer, X/X g 

X g thickness of the frozen layer at 

steady state 

x position coordinate in frozen layer 

measured from wall 

x’ dimensionless coordinate, x/X g 

y position coordinate in wall mea- 

sured from wall-coolant inter- 
face 

4 dimensionless coordinate, x/X 

p density of solidified material 

P w density of wall 

r time 

r' dimensionless time, 

Th i < T ! - T fV' >LX s 

Subscripts : 
c coolant 

f at freezing point 

ins insulated 

l liquid phase of solidifying sub- 

stance 
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w 


wall material 

test plate surface in contact with 
coolant 


2 


interface of solidifying material 
and wall 


1 


I,H, successive iterative approximate 
in solutions 
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APPENDIX B 


CORRECTION ON COOLANT TEMPERATURE TO ACCOUNT 
FOR HEAT LEAKAGE 


A cross section through the ice layer as shown in section B-B of figure 5(b) yields 
the configuration shown in figure 11. The heat flow q^ goes through the ice and the 
Inconel plate and corresponds to the one-dimensional heat flow considered in the theory. 
The heat flow qg is a heat leakage, and it is assumed to leak through the edges of the 
plate and does not pass through the plate in the manner of q^. This assumption was de- 
duced from some experiments in which the upper surface of the plate was insulated. For 
this condition, the temperatures T^ and Tg were found to be close together and were 

both above T m . 

c, m 

From equation (8), the heat flux at the ice-wall interface can be written as 


-PL — - h (T 
dr L 1 


Tf) + P c p 





(Bl) 


The heat flux from the plate to the coolant includes qg and is expressed by 

- + c *2 = h c( T l- T c ) m) 


t 2 -t. 


■w 


(B2) 


Equation (B2) is solved for T^ and substituted into equation (Bl). The result can then be 
placed in the form 


- p L^-VT*-T f >*Pc p | X ^ 


T -- IT c,m^ 


dx = - 




k w 


(B3) 


By comparing equation (B3) with equation (8), it is evident that the quantity T + (q 0 /h ) 
is the effective coolant temperature when heat leakage is present. For the tests where the 
water side of the test plate is insulated, Tg is equal to T^ and equation (B2) reduces to 
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‘>2' t >c( T l- Tc ,m) lns 

which provides an expression for qg. Thus, the effective coolant temperature that should 
be used when comparing the experiment with the one-dimensional analysis is 


T c 


= T 

c,m 



- T 

c,m 



(B4) 
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APPENDIX C 


CORRECTION OF COMPUTED HEAT FLUX TO COOLANT TO ALLOW FOR 
HEAT CAPACITY OF PLANE WALL 


For the wall on which the ice layer is growing (fig. 12), the heat-conduction equation 
is 


w 


3 2 T 3T 

w=p c —Z 

ay2 p w P,w 0T 


(Cl) 


Integrating twice with respect to y yields 


3T 


-k. 


w 


w 


3y 


= (T 2 - Tj) ± r r T dy dy 

a Z 1 a 3 t •'0 Jq w 


a r y 


y=o 


(C2) 


The term on the left is then expressed in terms of the convective heat-transfer coeffi- 
cient between the coolant and wall, 


k 


D C 


h c (T l ' T c ) =— (T 2 - Tj) - w p 


r 

3t -'o Jo 


T dy dy 
w J J 


(C 3) 


The double integral is the contribution from the specific heat term and is a correction for 
equation (45a). An approximate value for the double integral can be evaluated by using 
the following linear approximation for T w : 


T 2 ’ T i 

T w = T i + - 2 — 1 y 


(C4) 


After carrying out the integration, the heat flow to the coolant becomes 


h c (T l " T c ) = — ( T 2 “ T l> “ PwCp,Wa — ( 2T i + t 2 ) 


6 3t 


(C5) 
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The importance of the heat capacity correction can be evaluated from the measured T^(t) 
and TgCT) where the derivative 

i-[2 Tl (r) + T 2 (t)] 

is evaluated by finite differences. 
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Figure 1. - One-dimensional model of transient solidification in warm liquid 
flowing over convectively cooled plate. 
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Dimensionless frozen layer thickness, X/X, 


Dimensionless time, t‘ = xh z (Tj - T f )/pLX s 
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Figure 2. - Comparison of methods for predicting instantaneous thickness of frozen layer. (Some 
curves dashed for clarity only. ) 
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(c) Thermal resistance parameter, R « 1; subcooling parameter, S»3. 


Figure 3. - Concluded. 
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(a) Subcooling parameters, 0, 0.5, and 1. 

Figure 4. - Dimensionless growth times as function of thermal resistance parameter for various dimensionless thicknesses and 
subcooling parameters. 
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(b) Details of test section. 
Figure 5. - Test apparatus. 
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(a) Ice growth with initial ice patch nucleus on test plate. Thermal resistance parameter, R = 2.55; subcooling 
parameter, S = 0.232. 



20 40 60 80 100 200 400 600 800 1000 2000 4000 6000 8000 10 

Time from start of ice growth, t, sec 


(b) Ice growth with initial ice patch nucleus on test plate. Thermal resistance parameter, R - 1.67; subcooling 
parameter, S * 0. 225. 

Figure 7. - Comparison of measured with theoretically predicted ice thicknesses. (Plate temperature transients 
are in corresponding parts of fig. 8). 
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(c) Ice growth with initial ice patch nucleus on test plate. Thermal resistance parameter, R = 1. 35; 
subcooling parameter, S = 0.232. 
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(d) Ice growth with initial ice patch nucleus on test plate. Thermal resistance parameter, R = 0.631; 
subcooling parameter, $ = 0. 104. 


Figure 7. - Continued. 




(e) Ice growth without initial ice patch nucleus on test plate. Thermal resistance parameter, R = 1.21; 
subcooling parameter, S = 0.202. 
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(f) Ice growth with initial ice patch nucleus on test plate. Thermal resistance parameter, R = 1.29; 
subcooling parameter, S = 0.206. 


Figure 7. - Concluded. 



Temperature, 




Figure 8. - Temperature transients of test plate surfaces. (Each part of fig. 8 corresponds to same part of fig. 7. ) 
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(c) Thermal resistance parameter, R = 1.35; subcoofing parameter, S = 0.232. 
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(d) Thermal resistance parameter, R - 0.631; subcooling parameter, S - 0. 104. 
Figure 8. - Continued. 
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(e) Thermal resistance parameter, R = 1.21; subcooling parameter, S = 0.202. 
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(f) Thermal resistance parameter, R = 1.29; subcooling parameter, S = 0.206. 
Figure 8. - Concluded. 
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(a) Temperature-time history measured by probe at various fixed positions above plate surface. 
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(b) Temperature profile in water thermal boundary layer an instant before nucleation and 2 seconds later. 

Figure 9. - Temperature behavior in water thermal boundary layer before and after spontaneous ice nucleation on a subcooled plate. 
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